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Back-action ground state cooling of a micromechanical membrane via 

intensity-dependent interaction 
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We propose a theoretical scheme to show the possibility of achieving the quantum ground state 
cooling of a vibrating micromechanical membrane inside a high finesse optical cavity by back-action 
cooling approach. The scheme is based on an intensity-dependent coupling of the membrane to the 
intracavity radiation pressure field. We find the exact expression for the position and momentum 
variances of the membrane by solving the linearized quantum Langevin equations in the steady-state, 
conditioned by the Routh-Hurwitz criterion. We show that by varying the Lamb-Dicke parameter 
and the membrane's reflectivity one can effectively control the mean number of excitations of vibra- 
tion of the membrane and also cool down the system to micro-Kelvin temperatures. 

PACS numbers: 42.50.Lc, 42.79.Gn, 85.85. +j 
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I. INTRODUCTION 

The coupling of mechanical motion of a mechanical res- 
onator at the micro-and nano-meter scale to the electro- 
magnetic degrees of freedom via radiation pressure [lHa 
is widely employed for a large variety of applications [6[, 
more commonly as sensor to the detection of weak 
forces [3] and small displacements [8( or actuator in inte- 
grated electrical, optical, and opto-electronical systems^, 
[lOj . Modification of the resonator motion can be detected 
with high sensitivity by looking at the radiation (or elec- 
tric current) which interacted with the resonator. As an 
example, by measuring the frequency shift induced on 
the resonator one can detect the small masses. In prin- 
ciple, small displacements can be measured by detecting 
the corresponding phase shift of the light interacting with 
the resonator. 

Moreover, the rapid progress of nano-technology has 
enabled fabrication of micro- or nanomechanical res- 
onators [111 . Il2j with high frequency, low dissipation 
and small mass. In this direction the most experimen- 
tal and theoretical efforts [13l - [l9| are devoted to cool- 
ing and trapping such microresonators to their quantum 
ground state, but still in the attempt to observe quan- 
tized mechanical motion, the thermal fluctuation has be- 
come a major obstacle. The limited cooling efficiency 
and poor heat conduction at milli-Kclvin temperatures 
of cryogenic refrigerators have stimulated a number of 
studies on the active cooling of micro-mechanical res- 
onators in both classical and quantum regimes [g, l2Cf . 
Very recently, various types of optomechanical system 
have demonstrated significant cooling of the vibrational 
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mode of a mechanical resonator coupled to an optical 
cavity (171 l2ll |22| . Among them, the so-called membrane- 
in-the middle geometry |23|-[25| in which the mechanical 
degree of freedom is a flexible, partially transparent di- 
electric membrane placed inside a Fabry-Perot cavity 
with fixed end mirrors has attracted much attention. 
This has the advantage of not having to combine the 
flexibility needed for the mechanical oscillator with the 
rigidity of a high-finesse cavity mirror. Although the 
membrane is nearly transparent, it couples to the opti- 
cal cavity dispcrsively. This coupling is strong enough to 
laser-cool a 50-nm-thick dielectric membrane from room 
tcmpcraturc(294 K) down to 7mK|23j. 

From experimental point of view the cooling of a mem- 
brane has been achieved by exploiting in two different 
ways, back-action ground state cooli ng|15l Il7l. l25l . [26[ 
and cold-damping quantum feedback[2l|, [22|, [231 • By 
cold-damping quantum feedback, the oscillator position 
is measured through a phase-sensitive detection of the 
cavity output and the resulting photocurrent is used for a 
real-time correction of the dynamics, and with the back- 
action cooling, the off-resonant operation of the cavity 
results in a retarded back-action on the mechanical sys- 
tem and hence in a self-modification of its dynamics. In 
addition to this back-action effect, the radiation pressure 
couples mechanical oscillator to the cavity mode. Thus 
in appropriate conditions the later, plays as an effective 
additional reservoir for the oscillator. Therefore, strong 
radiation pressure coupling leads to significant cooling of 
the mechanical oscillator [18J. 

Motivated by the above-mentioned studies on the cool- 
ing of mechanical oscillator, in the present paper we 
deal with the study of cooling an optomechanical cav- 
ity with membrane-in-the middle geometry that consists 
of a high finesse cavity with two perfectly reflecting fixed 
end mirrors, and a partially reflective movable middle 
mirror(such as a dielectric membrane). In this type of op- 
tomechanical structure, the radiation pressure and cav- 



ity detuning arc periodic in the membrane displacement. 
This periodicity leads to an intensity-dependent inter- 
action between the intracavity radiation pressure and 
the membrane's motion. We show that the presence of 
intensity-dependent interaction modifies the effective me- 
chanical frequency and the effective mechanical damping 
of the membrane. It is also remarkable that in the steady- 
state condition, by controlling the system parameters, 
such as the membrane's reflectivity and the Lamb-Dickc 
parameter(LDP) the ground state cooling is approached. 
As we shall see by choosing proper values of these pa- 
rameters the system can be cooled down to micro-Kelvin 
temperatures. 

The paper is organized as follows. In Sec. II we de- 
rive an intensity-dependent Hamiltonian describing the 
coupling of a micromcchanical membrane to the radia- 
tion pressure field through j-phonon excitations of the 
vibrational sideband. In Sec. Ill we derive the quantum 
Langevin equations(QLEs) of the system and linearize 
them around the semiclassical steady-state. In Sec. IV 
we investigate the effective frequency and the effective 
damping parameter of the mechanical oscillator. In Sec. 
V the back-action ground state cooling is discussed. Fi- 
nally, a summary and some concluding remarks are given 
in Sec. VI. 



II. MODEL HAMILTONIAN 



fectly reflecting fixed end mirrors [23H25l |28|. In this type 
of optomechanical system the coupling between the mid- 
dle membrane and the optical cavity strongly depends on 
the position of the membrane. This position dependence 
results in a cavity detuning (for the fundamental mode 
of motion of the membrane and of the cavity) , which is a 
periodic function of the membrane displacement x, i.e., 
lu c (x) = (c/ L)cos~ 1 [\r c \cos( / Lirx/ A)] where L and r c are 
the cavity length and the field reflectivity of the mem- 
brane, respectively, and also the position of membrane is 
calculated from antinode of cavity field. The Hamilto- 
nian of the system is given by [2j, [24j 



H = tvjj c (x)a* a + hu> m b b + ihE(a < e 



iuilt 



" lt ), (1) 



where b and b^([b, tf] = 1) are the annihilation and cre- 
ation operators for the membrane with oscillation fre- 
quency u m and a and a^([a,a^] = l)are the annihilation 
and creation operators of the cavity mode with decay rate 
k. The last term in the above Hamiltonian describes the 
input driving by a laser with frequency uji , where E is re- 
lated to the input laser power Pq by \E\ = y / 2PoK/Huii. 
The fact that the motion of the membrane is quantized 
allows writing the operator of membrane's displacement 

x in terms of LDP, 



m = t 



(A and m denote, 

respectively, the wavelength of the incident field and the 
motional mass of the membrane) as 



We consider a high finesse cavity as is shown in Fig.l, x = T}o(tf + b). 

which is detuned by the motion of a partially reflective 
membrane placed between two macroscopic, rigid, per- By expanding tu c (x) in terms of b and tf we obtain 

I 



(2) 



(*) = ir " jl E E 



M r 



m=l fc=0 



(m-1)! rir, n (m-2k)(b+^) , h 1 

-irc m-i >~ l c i u.c.j-, 



(3) 



r 



where m = 21 + 1 and I is an integer number. By using without middle membrane, describes the free Hamilto- 
thc Baker-Campbell-Hausdorff theorem in Eq. (|3j) we may nian of the quantized cavity field and the free motion of 



rewrite the Hamiltonian of Eq.([T]) in the form 



the mechanical degree of freedom, and 



H = H + Hi nt + ihE(a^e' iuit - ae^), (4) H^ = hcJa[g*fj{n b )V + 9j (tfy 7> b )], (j = 0, 1, 2, ..) 



where 



Hq = hooker a + tko m b'b, 



(5) 



where 



9j — 2L \ l A V 2mu m / ' 



(G) 

(7) 



with wo = ttc/(2L) as the natural frequency of the cavity and 
I 



fj(n b ) = E E 



kc| r 



m=l k=0 



^ {e4faf (^^ S^gy ^ J(?)o)2(m _ 2fc)2]i 



(8) 



with n b = tfb and L? n as the associated Laguerre polyno- 
mial, describes a nonlinear coupling of the radiation pres- 
sure field with the movable membrane through j-phonon 
excitations of the vibrational sideband. 



III. LINEARIZATION OF QLES 

The dynamics of the system under consideration is gov- 
erned by the fluctuation-dissipation processes affecting 
both the optical and the mechanical modes. To study 
these effects, let us consider the first excitation of the vi- 
brational sideband by choosing j = 1 in the Hamiltonian 
(|4|) with redefining the motional annihilation operator 
b —> be™ and its conjugate 6^ — > b'e~~ m : 



H = fvuiod'a + huj m b'b — h ga' a[f (n b )b ■ 



-ihE{a ] t 



-iiUit 



tff{n b )] 



(9) 



where we have defined 

9 = 9\ 

f(n b ) : 



2ttc 
fi(n b 



2mu!„ 



(10) 



The nonlinearity function f(rib) determines the form of 
nonlincarity of the intensity-dependent of the coupling 



J 



between the cavity field and the membrane. We point 
out that in the limit of very small values of LDP, rjo 
and for certain values of the membrane reflectivity r c the 
nonlincarity function f(n b ) reduces to unity. Fig. 2(a) 
shows f(n b ) as a function of n b for 8770 and r c = 0.99, 
while Fig. 2(b) displays the behavior of f(n b ) versus n b for 
10 -4 77o and r c = 0.9, where we have set 770 = 10~ 5 . Ob- 
viously, for f(n b ) = 1 the Hamiltonian © reduces to the 
Hamiltonian of the standard opto- mechanical systemQ. 
Therefore, the inherent nonlinearity of the model under 
consideration can be attributed to the parameters rj and 



It should be noted that the experimental realization 
of the system under consideration shows 7?n << l,(e,g. 
for the experimental values given in Refs.|23M25l|. we ob- 
tain 770 < 10 -4 ). Therefore, one may keep terms up to 
first order in the phonons number n b in Eq.([5]) and ap- 
proximate the nonlinearity function ffj—i\(n b ) = f(n b ) 
by expanding the associated Laguerre polynomial: 



f(n b ) 



an bl 



(11) 



where we have defined the following real parameters 



c= E E (m-2k) 

m=l k=0 



JicE 



_. _ V V (»?o) 2 (»'-2fc) 3 \r c \" 
° ~ l^ 2^ 2! m 

ro=l fc=0 



771 iH^Dl ( e -i( I) o) 2 (m-2fe) 2 | 

k ) 4»-l[(2fi)I]^ C '' 



J 4™-M(^)!] at '• 



r 



(12) 



By substituting Eq. ([TTj) into the Hamiltonian © we 
obtain, 

H = H opt - hgaata(n b b + tfn b ) + ikE^e-^'* - ae iuit ), 

(13) 
where H opt = fiu)Q(Va + Huj m b'b — hgea)a{b + W) denotes 
the Hamiltonian of standard optomechanical system and 
the second term is the contribution associated with the 
dependence of the cavity detuning on the position of the 
membrane. As we see from the Hamiltonian (1 13|) the 
second term describes an intensity dependent coupling 
between the membrane's motion and the intracavity ra- 
diation pressure. 

Now, let us consider the fluctuation-dissipation theo- 
rem for the cavity field and the membrane's motion. For 
this purpose, We describe the effect of the fluctuations of 
the electromagnetic vacuum and the Brownian noise as- 
sociated with the coupling of the oscillating mirror to its 
thermal environment within the input-output formalism 



of quantum optics. For the given Hamiltonian (|13p in 
the interaction picture this yields the following nonlinear 
QLEs 



-(iAo + n)a + iga[e(b + b*) + <r(n b b + wn b )] 

+ E + V2Ka in (t), 



(14) 



b = -(iw m + 7)6 + igala[e + a(2n b + b 2 )} + y/2^b in (t), 

(15) 
where Ao = uj c — luq and 7 is the decay rate of the mo- 
tional phonons of the membrane. The cavity field quan- 
tum vacuum fluctuation ai n (t) satisfies the Markovian 
correlation functions 

(a m (t)al(t')) = [(n th ) + l]6(t-t'), 
(al(t)a in (t')) = (n th )6(t-t'), (16) 

(a ln (t)a m (t')) = (4 n (t)al(t')) = 0, 



with the average thermal photon number (nth)- In our 
study, we assume that the cavity is in zero temperature, 
i.e., (nth) = 0. It means that the number of thermal 
photons to be negligible at optical frequencies. Further- 
more, the motional quantum fluctuation bi n (t) satisfies 
the following relations 

(bin(t)bl n (t')) = [<n Mh ) + l]*(t-f), 

(b\ n (t)b m {t')) = (n b ,th)S(t-t'), (17) 

(hn(t)b in (t')) = (bt(t)bt(t'))=0, 

where (n^th) is the mean number of phonons in the 
absence of optomechanical coupling, determined by the 
temperature of the mechanical bath T 



(nb 



th 



1 



(18) 



Analyzing the quantum dynamics of the full nonlinear 
system is difficult, so we linearize the QLEs around the 
semiclassical fixed points. That is, we decompose each 
operator in Eqs. (fT4")) and (fl5|) as the sum of its steady- 
state value and a small fluctuation, e.g., b = b s + Sb and 
a = a s + Sa. This decouples our system into a set of 
nonlinear algebraic equations for the steady-state values 
and a set of QLEs for the fluctuation operators. We 
point out that the linearization of the QLEs is compatible 
with our aim i.e., the ground state cooling because as is 
shown in Ref.[3[, significant cooling is approached when 
the radiation pressure coupling is strong, which needs 
very intense intracavity fields i.e., |a s | >> 1. Thus the 
steady-state equations corresponding to the Eqs. (fT4")) and 
(fT5j) are given by 



x/k^+A 



t, bs = 



Ga s 



(19) 



where we have defined the effective detuning A and the 
effective coupling constant G as 

A = A -2gRe(b s )(e + a\b s \ 2 ) 1 (20) 

G = 2a s g[e + a{b 2 s +2\b s \ 2 )]. 

By eliminating the steady-state contribution and lineariz- 
ing the resulting equations for the fluctuations, we obtain 
the exact QLEs for the fluctuations of the quadrature op- 
erators 



5X 



A5Y - nSX 



2kX,. 



(21) 



SY = -ASX- K 5Y + (G R 5q + G I 5p) + V2KX m , 

(22) 



5q = fliSp-TiSq-dSX + y/2jq ir 



(23) 



S P = -n 2 Sq-T 2 6p + G R Sx + ^ Pm , (24) 



where Gr and Gi are, respectively, the real part and 
the imaginary part of G and we have symmetrized the 
fluctuation operators as 

8X=-=(5a + 5a s ),5Y = -—(6a -Sal), (25) 
V2 v 2« 

5q = -L(5b + Sb r ),6p= -^=(Sb-5tf). (26) 
V2 v 2i 

Furthermore we have defined 



* 'm 
fit 

r 2 



uj m - Agaa 2 s Re(b s ), 
n m + 2gaa 2 s Re(b s ), 

O m - 2gaa 2 s Re(b s ), 
7 + 2gcra 2 Im(b s ), 
7 — 2g<ra s Im(b s ). 



(27) 



We point out that in the limit of a — ► 0, Eqs. (|2Tjl - (f24j) 
reduce to the QLEs of standard optomechanical system 
[18| . These equations show that the nonlinear term in the 
Hamiltonian (fl3|) leads to additional parts in QLEs. In 
Eq.(H2), the term dSp and in Eq.flU]), the term -GiSX 
directly depend upon the parameter <r, which introduces 
further coupling in QLEs and this shows that the nonlin- 
earity is responsible for the appearance of the imaginary 
part of the effective coupling constant G. It is evident 
from Eqs. ([27|) that the modification due to the existence 
of nonlinearity is not limited just for G; the detuning and 
the damping parameters in Eg. ([27)) are also affected by 
the parameter a. 



IV. EFFECTIVE FREQUENCY AND 

EFFECTIVE DAMPING PARAMETER OF THE 

MEMBRANE 

In this section we evaluate the effective frequency oj e ff 
and the effective damping rate 7e// of the membrane in 
the system under consideration. For this purpose, we 
solve the linearized QLEs for the fluctuations in the dis- 
placement operator as 



where 



Mw) = xM/ T M, 



/tM^C^/tW, 



(28) 



(29) 



is the Fourier transformation of the fluctuations in the 
total force acting on the membrane, which includes a 
radiation vacuum and a Brownian motion component. 
Here x(w) describes the mechanical susceptibility of the 
membrane, given by 



X" 1 ^) = jf [ fi i°2 + (Ti - zcj)(F 2 + iuj) - niJ(w)], 



with 



I(u) 



_ A[G R +G I (r 1 -iLu)/n 1 ][G R -G I (r 2 -iLj)/n 1 } 



(fc-iw) 2 +A 2 -AGJ/Oi 



(30) 



(31) 



The mechanical susceptibility of the membrane can be 
considered as the susceptibility of an oscillator with ef- 
fective resonance frequency and effective damping rate, 
respectively, given by 



, ,2 (, .\ _ -n V J. O O A[^i+^ 2 ^ 2 -Gji^ 4 ] 

U ef} [bJ) - 1 il 2 + "l"2 - n 1 [K2 + (M-A') 2 ][K 2 + (cu+A') 2 ]' 

(32) 

7e//M .. _, , a[ M3 -(«+7)G?" 2 ] (-33) 

2 — I ~r JM/^ + ^-A'^H^ + tw+A') 2 ] ' 

where A' 2 = A 2 — AG 2 /Oi and we have defined 
Ml = ( K 2 + A' 2 )(r 1 G / + fi 1 G fl )(fi 1 G fl -r 2 G 7 ), 



m = (K 2 + r 1 r 2 -2K{r 1 + r 2 ) + A ,2 )Gj 

+ {T 2 -T 1 )n 1 G I G R -nlG 2 R , 

Ma = 7 G 2 ( K 2 + A' 2 ) + 

fi 2 G| + n 1 G I G R (T 1 - r 2 ) - G^r* 



(34) 



The modification of the mechanical frequency due to the 
radiation pressure shown by Eq. ([32)) is the so-called opti- 
cal spring effect |18| . The existence of nonlinearity in the 
Hamiltonian (fT3|) leads to the appearance of higher order 
of system response frequency u>, in contrast to the re- 



sults of the linear coupling theory [18(- It is evident that 
in the limit of u — > the terms Gjuj 4 /rt 2 in Eq. (|32]) 
and (k + j)(Gj/Qi) 2 uj 2 in Eq. ([33|) are removed. Under 
this condition the effective frequency and the effective 
damping parameter reduce to the corresponding parame- 
ters in standard optomechanical system(See Eqs.(18) and 
(19) in Ref . fig ) . Fig. 3(a) shows the normalized effec- 
tive frequency as a function of the normalized system 
response frequency io/uj m at A = u> m for different values 
of LDP. As is seen, LDP alters significantly the effective 
frequency. The normalized effective damping rate ver- 
sus the normalized system response frequency w/w„, at 
A = w m have been plotted in Fig. 3(b). As is seen, by in- 
creasing LDP the effective damping rate increases. The 
significate increasing of the effective mechanical damp- 
ing rate is the basis of the cooling process. The ground 
state cooling can be obtained only if the initial mean 



thermal excitation number n = [exp(huj m / ksT) — 1] _1 is 
not prohibitively large. This can be approached for large 
values of the effective mechanical damping rate. As is 
clear from Fig. 3(b) for positive A and by choosing proper 
values of LDP the effective mechanical damping rate is 
significantly increased. 



V. BACK-ACTION GROUND STATE COOLING 

IN THE PRESENCE OF 

INTENSITY-DEPENDENT COUPLING 

In this section we study the effective ground state cool- 
ing of the membrane's motion. We first evaluate the 
mean energy in the steady-state by defining the system's 
potential as 



U 



^((6q(t) 2 ) + (Sp(t) 2 )) = fi Wm (n e// + \). (35) 



The ground state is approached for n e ff ~ or U — 
fko m /2. This is obtained if (5q(t) 2 ) ~ (Sp(t) 2 ) ~ 1/2 
occurring in steady-state. The system reaches a steady- 
state only if it is stable. In this situation all the poles 
of the effective susceptibility x lie i n the lower complex 
half-plane. The stability conditions can be attained by 
applying the Routh-Hurwitz criterion [29( which gives the 
following stability expressions 



Mi 



M 2 



( K 2 + A 2 )(n 1 n 2 + r 1 r 2 ) 

-A[n 2 GJ + OiG| + GnGjir, - r 2 )] > 0, 



(36) 



A 4 + A 2 [ 7 2 - 2L 1 L 2 + 4k 7 + 2{k 2 - fiifi 2 )] 
2(k + 7 ) 2 



+A- 



IV) 



-p 2 Gj + n x G R + G R G I {T 1 -V 2 )] 



+[nin 2 + (Ti + «)(T 2 + n)] 2 > 0. 

Now, we turn our attention to Eq. (|28p in which /p(w) 
is the fluctuations in the total force acting on the mem- 
brane's motion and it has the following form 



J 



/ T (w) = V2 7 p in + 






'2n 



r 2 - iu - 

QiGr — 



ApxGn - Gj(T 2 - iu)] 



CIi(k - itu) 2 + iliA 2 
Gi{T 2 -ilo) 



f2i(/-c 



fiiA 2 - AG 2 



AG 2 . 
AY in + (k - iu>)Xi. 



(37) 



r 



where the first and second terms show the fluctuations 
in the membrane's motion, and the third term describes 
the fluctuations in the intracavity field. By using Eq.(| 



we obtain 



(6q(u)Sq(u>')) = xMxO^X/rM/iV)), 



(38) 



and by inverse Fourier transforming Eq. (138|) we have 
1 



(Sq(t)Sq(t>)) = 



2tt 



oo /»00 



dbJ I did X 

— oo J —oo 

— iu)t Auj't' 



(39) 



e -«, * e «-* x ( w ) x ( w ')(/ r («)/r(w'))- 



By using Eq. (|37p and the quantum fluctuation relations 
([TB| and (HH we find 



(J T (u)f T (u)')) = j(2n b , th + l) + 



l{2n b ,th + 1) 



+ 2-f(u + AIm[ 



CM 



fii 



r 2 - iuj - 



AC(u) 



C(u) 



D(w) 



D(w) 



(A 2 + k 2 + w 2 + Aw) <y(w + w'), 



(40) 



r 



where we have defined 

C(u) =G R -G I (T 2 -iu)/Q 1 , 

D{uj) = ( K -iu) 2 + A' 2 . (41) 

The integral of Eq. ([39| for the displacement variance can 
be solved exactly when the stability conditions given by 
Eqs. (|36p are satisfied. Performing the integral and using 
x(— w) = X*( w ) w ith setting t = t' , one gets the final 
expression for the time dependent position variance 



(Sq 2 ) = j(2n b ,th + l)JVi + kN 2 



(42) 



where 



2N X = n 



2 (a 1 b 2 -a 3 b 1 b 3 (a 3 - 01O2) 



K 



M 2 M 2 Mt 

aid 2 - a 3 di + a 2 a 3 - aia 4 ^3(03 — O1O2) 



27Vo 



O 



M 2 
2/ Q i c 2 - Q3C1 c 3 (a 3 — OiOg) 
M 2 M 2 Mi 



M 2 Mi 



;(' 



(43) 



By applying the same procedure, we obtain the following 
expression for the momentum variance 



(Sp 2 ) = 7(2n t)tfc + 1)N[ + kN 2 , 



(44) 



where 



2N{ = fi 2 ( 



Oi6 2 - a 3 bi 63(03 - aia 2 ) 



M 2 M 2 Mi 

aid 2 — 03d! + a 2 a3 — 0104 d 3 (a 3 — a\a 2 ) 



2N 2 ---- fi 2 ( 



M 2 
, ( aic 2 - a 3 c[ c' 3 {a 3 - aia 2 ) 
M 2 M 2 M X 



M 2 Mi 



(45) 



The explicit expressions for the parameters ai,bi,Ci, di 
and b'i, c^, c^ are given in the appendix. 

The effective phonon number of the membrane's mo- 
tion can be defined from Eq. ([35| 



n e// = \{{5q{t?) + (Sp(t) 2 



1) 



As we mentioned, the ground state cooling is approached 
if H e // < 1 which is not the only condition for cool- 
ing. In order to get the quantum ground state cool- 
ing it is also necessary that the uncertainly principle 
is stratified for both momentum and displacement vari- 
ances. This means that both variances have to tend 
to (Sq 2 ) ~ (Sp 2 ) — 1/2, therefore energy equipartition 
has to be satisfied in the optimal regime close to the 
ground state. We point out that, in general, one does 
not have energy equipartition because these variances are 
not equal (5q(t) 2 ) ^ (Sp{t) 2 ). 

In order to have an intuitive picture, we have plotted 
the both variances in Fig. 4 for different values of LDP 
and for the experimentally feasible parameters (see Ta- 
bid in Ref.[23| and also [2J| ), i.e., a cantilever with 
motional mass m = 50pg,u> m /2ir = 10 5 Hz, r c = 0.999 
and Q — 1.2 x 10 placed inside an optical cavity with 
length L = 6.7cm and damping rate k = 0.047w Tn driven 
by a laser with wavelength Ao = 1064nm and power 
P c = 50^1^. It is evident from Fig. 4 that by varying 
LDP one can access the ground state cooling condition 
(Sq 2 ) ~ (Sp 2 ) ~ 1/2. As is seen from fig.4(c)for ?7o around 
A = l.lWm one has (Sq 2 ) ~ (Sp 2 ) ~ 1/2. This means 
that LDP may significantly alter the cooling process and 
it can be used as an additional control parameter to op- 
timize the best condition for the ground state cooling in 
the system. 

Fig. 5(a) shows the effective phonon number as a func- 
tion of detuning for three different values of LDP. It is 
clear that with the increasing value of LDP the mini- 
mum value of n e ff is shifted toward zero, such that for 
770 w e have Ji e // — 0- To investigate the effect of the 
membrane's reflectivity on ground state cooling we have 
plotted the effective phonon number as a function of de- 
tuning in Fig. 5(b). It is interesting that increasing of r c 
significantly decreases the effective phonon number. One 
can also define an effective temperature with respect to 
the effective mean excitation number n e // as 



T, 



huJr, 



e ff k B Zn(l + l/n e// ) " 



(47) 



(46) we have plotted the effective temperature for experimen- 



tal values given in Fig. 4 and for the initial reservoir tem- 
perature T = OAK corresponding to n ~ 83306 against 
the normalized detuning and for different values of LDP. 
It is remarkable that the minimum achievable tempera- 
ture is T e ff ~ \[iK . Therefore by controlling LDP and 
r c one can control the minimum value of the effective 
temperature of the system. 



VI. SUMMARY AND CONCLUSIONS 

In this paper we have introduced a physical scheme to 
develop the ground state cooling of a mechanical mem- 
brane. As we have seen, the dependence of cavity detun- 
ing to the membrane displacement leads to an intensity- 
dependent interaction between the intracavity radiation 
pressure and the membrane's motion. Consequently, an 
additional nonlinear term is appeared in the Hamiltonian 
of standard optomechanical system. Existence of this 
part in the Hamiltonian alters the dynamical behaviour 
of the system by introducing extra coupling terms in 
QLEs. We have derived the effective mechanical suscep- 
tibility and have studied the effect of LDP on the effective 
mechanical frequency and the effective damping param- 
eter. 

As we have also seen, the additional nonlinear term en- 
hances significantly the back-action ground state cool- 
ing process to obtain smaller effective temperature. By 
varying LDP and membrane's reflectivity one can op- 
timize the cooling process to achieve the ground state. 
We have found that for initial reservoir temperature 
T = 0.4?f (corresponding to n ~ 83306) and in the exper- 
imentally accessible parameter regimes, one can observe 
the mean mechanical excitation number close to zero 
with corresponding effective temperature around micro- 
Kelvin. 
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Appendix: solution of rational integral of Eq. (|39p 

when the stability condition is satisfied, according the 
Routh-Hurwitz criterion one can solve exactly the inte- 
gral of rational function of Eq.(|39|)to find the position 
and momentum variances of the membrane given, respec- 
tively, by Eq.lO and by Ea.([4"4"|) in which 



<"l 



2( 7 + K ),a 2 = -(A 2 



T 1 r 2 + n 1 n 2 + ^K), 



= -2[7(A 2 



K (r 1 r 2 + r! 1 r! 2 )], 



1, 6 2 = 2(k 2 -A 2 ), b, = (K 2 + A 2 3 ) 2 , 

,2 i 2\„ , (n _ GlT 2s2 



(^L) 2 ,c 2 = (A 2 + K 2 )c 1 + (G R 



Ol 



c 3 = (A 2 + K 2 )(G 



G I^2 ]2 



(A.l) 



di 

d 2 ■- 

da = 



and 



2(k 2 - A 2 ) + r 2 . 

(fee 2 + A 2 ) 2 + 2T 2 2 (kc 2 - A 2 ) + 2(2fcc + r 2 )AG/G fl , 

(k 2 + A 2 )r 2 - Add 



b' 2 = 2{n 2 -A' 2 ), b', = {n 2 + A'if, 

A = (g) 2 , C2 = (A 2 + K 2 K + (G / + ^i) 2 , 



(A 2 + K 2 )(d + 



GrTi 



d' 2 = 
d' = 



n 2 

J2 i a 2\ 2 i ot^2 / ;„„2 



(A.2) 



(kc z + A Z Y + 2T{{kc z - A 2 ) - 2(2fcc + Ti)AdG R , 

( K 2 + A 2 )r 1 + AG / G 1 1 ' 



with A' 2 = A 2 - 4^ 
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Figure Captions 



Fig.l: A high- finesse optical cavity with two rigid 
end mirrors and a dielectric membrane centered at the 
antinode of cavity field. 



/(?!{,) as a func- 
= 0.99, 8?7o- (b) 



Fig. 2: The nonlinearity function 

tion of phonon number rib for: (a) i ^ 

r c = 0.9, 10~ 4 ?7o. Here we have set rj — 10 . 

Fig.3:(a)The effective mechanical frequency of Eq. ([52"l) 
and (b) The effective mechanical damping rate of 
Eq. (I33p versus normalized frequency. Parameter values 
are!!!!! W m /27r = W 5 Hz,L = 7cm, P = 61 fiW, K = 
0.051cj m ,Q = 1.2 x 10 7 , motional mass m = 0.5pg and 
770 = 6.8 x 10~ 7 (i)0.65?7o(red line),(M)0.7?7o(green line) 
and (m)rjo(blue line). 

Fig. 4: the membrane's position variance (<5g 2 )(red line) 
and the membrane's momentum variance (Sp 2 ) (blue 
line) versus the normalized effective detuning for the 
initial reservoir temperature T = 0.4 corresponding to 
n e ff ~ 83306 and for different values of LDP:(a)0.85?7o 
(b)0.977o and (0)770. 

Fig. 5: The effective mean excitation n e ff versus 
the normalized effective detuning for the initial reservoir 
temperature T = 0.4 corresponding to n e // — 83306 
(a)for different values of LDP: 0.85?7o(red line), O.9?7o 
(blue line) and 770 (green line)(b) for different values of 
membrane's reflectivity: r c = 0.98(red line), r c = 0.99 
(blue line) and r c = 0.9999(grcen line). 



Fig. 6: The effective temperature of Eq. pTj) versus 
the normalized effective detuning for different values of 
LDP: O.85?7o(red line), O.9?7o (blue line) and r7 (green 
line) . 
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